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1. Introduction 

Loop models are a general class of statistical problems where the elementary statistical 
objects are loops drawn on a two-dimensional lattice (for a review see for instance [0). 
Loop models arise naturally in the high temperature expansion of lattice statistical models 
but also as the description of one-dimensional lattice objects, like self-avoiding polygons. 
An interesting limiting case in that of fully packed loop models where the set of loops is 
required to cover the entire lattice, without vacancies. By assigning an activity n per loop, 
such models can be thought of as the zero temperature limit of 0{n) models, more simply 
referred to as fully packed 0{n) models. In the following, we shall restrict our discussion 
to triangular lattices, be it the flat regular triangular lattice or random triangulations. 
The loops are understood here as being self-avoiding and non-intersecting. Furthermore, 
all the triangles of the lattice are assumed to be visited by a loop. For n = 1, the fully 
packed 0{n) model describes for instance configurations of dimers drawn on the dual of the 
(regular or random) lattice, the dimers occupying the edges dual to those not traversed 
by a loop (each triangle has exactly one such edge). The model describes equivalently 
the ground states of an anti-ferromagnetic Ising model. In the limit n ^ 0, the model 
describes Hamiltonian cycles on the lattice, which are the compact conformations of a 
polymer ring. At n = 2, the model defined on the regular triangular lattice is equivalent 
to a three-coloring problem, namely the problem of coloring the edges of the lattice with 
three colors in such a way that the three edges adjacent to any triangle are of different 
color 0. Alternatively, it describes the possible folded states of the regular triangular 
lattice onto itself 0. On random triangulations, this equivalence with three-coloring and 
folding problems remains valid only for a restricted class of triangulations [^^ . We shall 
return to this below. 

A remarkable prediction for the fully packed 0{n) model on the regular two- 
dimensional triangular lattice is that it is not in the same universality class as the usual low 
temperature dense phase of the 0{n) model in which vacancies are allowed. If we denote 
by Cd{n) the central charge of the dense phase fixed point, that of the fully packed phase 
is expected to present a shift by one, namely c/(n) = Cd{n) + 1. This remarkable fact was 
first conjectured in on the basis of transfer matrix studies, and then confirmed in [^ on 
the grounds of a nested Bethe Ansatz solution. As explained in [Q , the shift by one in the 
central charge can be given a nice heuristic interpretation. Indeed, by marking those edges 
of the triangular lattice which are traversed by the loops, any set of fully packed loops 
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can be viewed as a two-dimensional picture of a three-dimensional piling of cubes, whose 
surface defines a one-dimensional SOS height variable on the triangular lattice. This SOS 
degree of freedom, which emerges only if the loops are fully packed, is responsible for the 
shift in the central charge. 

Unfortunately, the above geometrical picture breaks down when going to ordinary 
random triangulations. Indeed, the local rules which would define a height variable out of 
the loops in general lead to frustrations. In this case, the SOS variable cannot be properly 
defined anymore and the shift in the central charge does not occur. In other words, for 
ordinary random triangulations, the fully packed 0{n) model is again described by the 
dense phase fixed point. 

Recently, in [Q, it was conjectured that the shift by one in the central charge can 
be reinstated in the random case if the triangulations are restricted to the class of so- 
called Eulerian triangulations. An Eulerian triangulation is a triangulation where an even 
number of triangles meet at any given vertex. Eulerian triangulations arise naturally in 
the context of folding problems involving random lattices. Indeed, in genus zero, Eulerian 
triangulations are the vertex-tricolorable triangulations, namely those for which each vertex 
can be assigned one of three colors in such a way that any two neighbors have distinct 
colors 0. This latter condition ensures the possibility of folding the triangulation in 
two dimensions, as explained in Q. For Eulerian triangulations, the different possible 
folded states can then be mapped onto edge-three-colored states, or equivalently onto 
configurations of fully packed loops with a weight n = 2 per loop. The fully packed 
0{n = 2) model on random Eulerian triangulations thus provides a natural description 
of the folding of fiuid membranes. As explained in H], for Eulerian triangulations, the 
construction of the SOS variable again becomes possible without frustrations, and a shift 
by one in the central charge should then be observed. This phenomenon was moreover 
confirmed numerically in [§ in the limit n — > by a direct counting of Hamiltonian 
cycles on random Eulerian triangulations with up to 40 triangles. The string susceptibility 
exponent was found to be compatible with the value 7 = (— 1 — vT3)/6 expected for a 
central charge c = — 1, instead of the value 7 = — 1 found for ordinary triangulations, 
corresponding to Q(n = 0) = — 2 []T^-|T^|. 



The purpose of this paper is to confirm this shift phenomenon in the case n = 1 by 
showing that the fully packed 0{n = 1) model has c = 1 when defined on random Eulerian 
triangulations, as opposed to the usual result c = for the (dense or fully packed) 0(n = 1) 



model on ordinary random triangulations [13,14 



The paper is organized as follows: in Sect. 2, we present a matrix model formulation 
for the 0{n = 1) model on random Eulerian triangulations. This model is shown in Sect. 3 
to be equivalent to a particular gravitational 6-vertex model, described at criticality by a 
c = 1 conformal field theory coupled to gravity. In Sect. 4, we extend our matrix model 
formulation to arbitrary values of n. We focus in particular on the limit ti ^ describing 
Hamiltonian cycles. A few concluding remarks are gathered in Sect. 5. 

2. Matrix Model for the Fully Packed 0{n = 1) Model on Random Eulerian 
triangulations 

As mentioned in the introduction, an Eulerian triangulation is a closed random trian- 
gulation of arbitrary genus, for which each vertex has an even number of adjacent triangles. 
Alternatively, an Eulerian triangulation can be defined as a triangulation where one may 
associate a sign + or — to each triangle in such a way that any two adjacent triangles have 
opposite signs. 

Here we consider random Eulerian triangulations equipped with fully packed self- 
avoiding loops of adjacent triangles, i.e., triangulations covered by a set of loops such that 
each triangle belongs to exactly one loop. In this section and in Sect. 3, we address the 
case of an activity n = 1 per loop; the case of general n (including the Hamiltonian cycle 
limit 71 ^ 0) will be discussed in Sect. 4. As opposed to the usual 0{n) model coupled to 
gravity [|l^, we insist on imposing here the two crucial restrictions: (1) the loops must 
be fully packed (no vacancies) and (2) the triangulations must be Eulerian. Only in this 
particular case is the model expected to be described by a different universality class than 
the usual 0{n) model, with the c -^ c + 1 shift in the central charge. 

Our model is best expressed in the dual picture as that of a three-coordinate {(/}^) 
lattice with bi-colored vertices (corresponding to the above-mentioned -|- and — signs), 
equipped with loops visiting all vertices. In particular, the signs of the vertices visited 
by a given loop alternate along the loop. The corresponding graphs are the Feynman 
diagrams of the following simple matrix model. We consider a pair of complex N x N 
matrices (X, L) , where (X, L) will correspond to the -|- vertices of the dual graph, whereas 
(X^jL^) will correspond to the — vertices. As usual, Feynman diagrams for such objects 
are obtained by joining pairs of double-edges, each pair corresponding to a matrix element 



Mij (resp. MM, where the two hnes carry the matrix indices i and j, and the double-edge 
is oriented away from (resp. towards) a vertex. We need the following interactions 

X 



Tr(XL^) : 


L L 




X' 


Tr(Xt(Lt)2) : 


i e 


and propagators 





(2.1) 



X X^ 

{{X)ij{X^)ki) = ^Sii5jk = ]:::::::):::::::[ 

L L^ 

(2.2) 
The latter form the edges of the ^^-diagram, and are of two types: those visited by the 
loops {LL^), and those not visited by loops {XX^). The mixed nature (MM^) of the 
propagators guarantees the Eulerian structure. The interaction terms ( |2.1|) describe the 
three-coordinate vertices with the corresponding signs and exactly two occupied (L) and 
one empty {X) incoming edges. The diagrams with vertices ( |2.1| ) and propagators ( |2.2| ) 
arise in the Feynman expansion of the following four-matrix integral 

y (2.3) 

V{X,L) = XX^ + W - g{XL^ + X^{L^f), 

where the standard measure over N x N complex matrices reads dMdM^ oc Y[i<i j<n 
dKe{Mij) dlni{Mij), and is normalized in such a way that Z{g = 0; A^) = 1. 

As usual, the free energy f{g]N) = Log Z{g]N) = Yl,h>o^'^~'^^fh{9) is expressed as 
a sum over the contributions of the connected Eulerian triangulations of genus h. The 
genus zero limit is therefore obtained by taking N ^ oo. We have represented in fig.l an 
example of a connected genus zero diagram with eight triangles and two loops. 




Fig. 1: A typical genus zero configuration involving two loops (solid double- 
lines) fully packed on an Eulerian triangulation made of 8 triangles. We have 
represented by dashed double-lines the un-occupied edges of the dual lattice. 
The orientation of the double-lines reflects the Eulerian constraint (all arrows 
point towards triangles with -t- signs, and away from those with — signs, hence 
the orientation alternates along each loop). 



3. Mapping to a critical point of the gravitational 6-vertex model 



The integral ( |2.3| ) is Gaussian in all matrices. Let us first integrate over X by setting 
X = -y=(P + iQ)j where P and Q are two N x N Hermitian matrices, and the measure is 
transformed into dXdX"^ oc dPdQ, where dP and dQ stand for the standard Haar measure 
for Hermitian matrices, normalized in such a way that Z{g = 0; A^) = 1. Similarly, we set 
L = -y=(A + iB), with A and B Hermitian, so that the potential becomes 

Tr{V{X,L)) = TrQ(A2 + s2) + l(p2^g2^ 

- -^ (P(A2 - B^) - Q{AB + BA))] 
v2 / 

= Tr(^{A^ + B^) (3.1) 

+ l{p-j={A^-B^)r-'-^iA^-Br 

+ i(Q + -^{AB + BA)f - ^{AB + BAf 



2'^ y/2' " 4 



Performing the Gaussian integrals over the shifted matrices P and Q, we are left with 



Z{g;N) 
W{A,B) 



rfArfSe-^T--^^^'^), 



1 



{A^ + S^ 



^^.4 



(3.2) 



(A^ + 5^ 



iABf 



This is nothing but the partition function of the gravitational 6-vertex model solved in 
the large- A^ limit by Kazakov and Zinn-Justin [16|| (with the parameters a = P = g'^), 
whose critical point g = Qc = l/(2-\/7i") corresponds to a compactified boson with radius 
R = 1/{2V2). 

The crucial outcome of this equivalence is that the conformal field theory underlying 
our problem has central charge c = 1, which precisely corresponds to a shift by one 
of the central charge c = of the ordinary (fully packed or not) dense phase of the 
0{n = 1) model on arbitrary (i.e. non-necessarily Eulerian) triangulations. This proves in 
the particular case n = 1 the claim of [^ that the central charge increases by one when the 
fully packed model is defined on random Eulerian triangulations as opposed to ordinary 
random triangulations. 

Two remarks are in order. From the critical value gc — l/{2^/n), we deduce that the 
number of genus zero Eulerian triangulations with 2T triangles and equipped with fully 
packed loops behaves for large T as 



.E 



Z2T{n = 1) 



(47r) 



(3.3) 



to be compared with z^ ~ 8^ for pure Eulerian triangulations 0, to -2^(n = 0) ~ 
(10.10...)"^ for Eulerian triangulations equipped with Hamiltonian cycles f^, and finally to 
Z2T{n = 1) ~ (24)^ for ordinary random triangulations equipped with fully packed loops. 
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Fig. 2: Shrinking the (XX^) propagators (dashed double-edges) produces a 
particular 4-valent vertex of the 6-vertex model. 
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Secondly, let us note that the 6-vertex correspondence is best seen by shrinking the 
(XX^) propagators so as to form four-valent vertices with oriented edges (cf. fig. 2). 



4. The 0{n) model on random Eulerian triangulations 

4.1. Matrix model for arbitrary n 

The fully packed 0{n) model must incorporate a weight n per loop, obtained for 
integer n by replicating n times the N x N matrix L of ( p.3|) . Beside the complex matrix 
X , we therefore introduce n complex matrices L^, a = 1, 2, ..., n, with the following vertex 
interactions 

Tv{XLl) and Tt{X^{LI)^), for a = 1, 2, ...,n, (4.1) 



and propagators 



{{X)ij{X^)ki) = —SiiSjk, 

^^ (4.2) 



■t' 



1 



{{La)ij{Lp)kl) = —Sa(3Sil5jk. 

This allows only La-matrices of the same color a to form loops. The corresponding matrix 
model partition function reads 



-iVTr(l/(X,Li,...,L„)) 
a = l 



/ll. 
dXdX^ Y[ dL^dLl e" 

a = l 
n 

y(X,Li,...,L^) = XXt + ^L^Lj^ (4.3) 

a=l 
n n 

-g{xj2iL^r+x^j2(Lir)- 

Here again, the integration measure is normalized in such a way that Z{n, g = 0; N) = 1, 
and the net result in the perturbative expansion of Z{n, g; N) is to attach a weight n per 
loop of L-matrices. 

Contrary to the n = 1 case, let us first integrate over the n matrices Lq. = (Aq, + 
i-BQ,)/v2, where Aa and Sq,, ct = 1, 2, ...,71 are n Hermitian matrices of size N x N. To 
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do this integration, we note that the potential V of (^4.3| ) takes the form 

ct 

- |X ^ {Al ~Bl+ i{A^B^ + B^A^)) 

a 

- |xt Y, {Al -Bl- z{A^B^ + B^A^)) 



(4.4) 



a=l ^ ^ 



where the quadratic form Q reads 
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Q = {I I)l2 - ^{{X I + I X)K + {X^ I + I X^)K). (4.5) 

Here we have denoted by / (resp. I2) the N x N (resp. 2x2) identity matrix, and K, K 
are the foUowing 2x2 matrices 

Performing the Gaussian integration over the A's and S's, we finaUy get 

Z(n,(7;iV) = /"rfXc^Xtdet(Q)-"/2e-^^'-(^^'^ (4.7) 

We may further expand 

det(Q)-"/2 = exp[-^Tr Log Q] 

°° 2m (4.8) 

= exp[n^ ^—-Tx[{{X®I + I®X)[X^®I + I®X^)Y")\ 

where we have used the fact that K^ = K^ = 0, hence only the terms of the form 
Tr[(i^K)"'] = 2^"^ or Tr[(i?ii:)"'] = 2^"^ contribute, which cancel the 1/22"^+^ pre-factor. 
This reduces the 0{'n) model partition function on Eulerian triangulations to a Gaussian 
complex one-matrix model with some specific integrand. Note that in eqns. ( [4.7| ) and ( [4.8| ) 
the parameter n can now take any real value. 



4-2. The n ^ limit: Gaussian matrix model 

The result ( [4.7|) yields in particular in the limit n — » the partition function for 
Hamiltonian cycles on Eulerian triangulations, expressed as a complex one-matrix integral 

= 5^|^(Tr([(X®/ + /®X)(Xt®/ + /®Xt)]'")) 

7Tl=l 

OO 2m 

= y ^— V {Tr(X''UX^Y^X''\.m(X^-''UX^V-''^X^-''^..)), 

m=l !/i,...,!y2™e{0,l} 

(4.9) 
where the bracket stands for the Gaussian integration over the complex matrix X, namely 
(/(X)) = JdXdX^f{X) exp{-NTT{XX^)), and is normalized in such a way that (1) = 1. 
Moreover, the large- A^ limit of ([4.9|), znig) = hniTv^oo -^Zuig^N) yields the generating 
function for genus zero Eulerian triangulations equipped with Hamiltonian cycles. Due to 
the known large- A^ factorization property (Tr(/(X)^(X))) ~ (Tr(/(X)))(Tr(^(X))), we 
also have 

oo 2m 

^^(^) = E|^ E {{Tr{X^^{X^Y'X-K.))){{TT{X^-^-{X^f-^-X^--K.))), 

m=l i/i,...,!y2^6{0,l} 

(4.10) 
where the double bracket is defined as ((Tr/(X))) = limAr^oo ;^(Tr/(-^))- 

The formula ( |4.9| ) can be interpreted pictorially as follows. The quantity \{X ® I + 
I ® X){X'^ ® I -\- 1 ® -^^)] can be represented as a succession of 2m points along a line, 
from which oriented bonds originate, alternately oriented away from and towards the line 
to account for the alternation of X and X^, and going above (resp. below) the line if a 
term M ® I (resp. / ® M), M = X, X^ is selected i.e. according to whether Vi = 1 (resp. 
Ui = 0) at the z-th point in ( 4.10| ). This gives rise to simple diagrams like that of fig. 3, 



where the 2m points are connected pairwise by oriented arches either above or below the 
line, such that arrows inwards and outwards alternate along the line. This alternation can 
be replaced by a sequence of alternating + and — signs, each arch connecting a -f to a — 
as in g. 

In ( [4.10| ) , we see that the knowledge of all Gaussian averages of traces of words in X 
and X'^ would immediately give access to znig)- Such objects have been considered in |1T5[] , 
as Gaussian averages of traces of words involving Hermitian matrices, for which a complete 

9 



/\ \/ 



/\ 



XX^<^I 



\/ /\ 



/\ 



x^x^ 



\/ 



\/ 




I^XX^ x^<^x 

(a) 



(b) 



Fig. 3: (a): rules for representing a term in the expansion (|4.10|) . The four 
cases correspond respectively to the values (1^1,2^2) = (H), (10), (00), (01), 
namely to the selection of X^i(Xt)^2 X'^-''^{X^y-''^ in the product {X 
/ + / (g) X){X^ / + / ® X^). It is understood that the line representing 
X should be drawn to the left of the line representing X^ . In the figure 
(b), we have represented a typical term, with 2m = 14 points, and a choice 
(z^i...z^i4) = (10011110011001). The value u = 1 (resp. 0) corresponds to an 
arch going above (resp. below) the line. The Eulerian condition imposes that 
orientations alternate between successive points. 



set of recursion relations has been found, solving in principle (but unfortunately not in 
practice) our problem. These were studied in the context of meander enumeration, namely 
the enumeration of possibly interlocking loops (roads) crossing a line (river) through 2m 
given points (bridges). The generating function for meanders with a weight 2 per connected 
component of the road can actually be recast in a way very similar to (|4.9|) , namely 



M{g^;N) 



771=1 



Jlra 
2171 



Y^ |— (Tr((X®Xt+Xt®X)'")). 



(4.11) 



This corresponds precisely to retaining in (|4.9| ) only the sets of z^'s that satisfy i'2i = 
l — U2i-i for i = 1, 2, ..., TO which again corresponds to retaining only the two configurations 
on the right hand side of fig. 3- (a). These two can both be recombined into a single oriented 
bond crossing the line (with a weight g^ per intersection), leading to the usual picture 
of a multi-component meander with two possible orientations per connected component. 
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accounting for the factor of 2. In the planar (A^ -^ oo) hmit, this yields 

oo 4?Tl 

^(^') = Ei;^ E ((Tr(X-(Xt)i--X-(Xt)i--..))) 



X ((Tr(X^-^i(Xt)'^iX^-'^2(Xt)^2_.)))- 

(4.12) 



5. Conclusion 

In this paper, we have considered some fully packed loop models on Eulerian tri- 
angulations. In the case of the 0{n = 1) model, we have shown that taking Eulerian 
triangulations rather than arbitrary ones leads to a shift c ^ c + 1 in the central charge of 
the conformal theory coupled to gravity describing the corresponding critical point. More 
generally, given any matrix model describing random triangulations, typically defined by 
a Hermitian multi-matrix integral with a potential of the form 



1 ^ 

y(Ai,...,Ap) = -J2A^^-J2c,JkA^A^Ak, (5.1) 



i=l ijk 

we can restrict ourselves to the class of Eulerian triangulations by replacing the Hermitian 
matrices Ai by complex matrices Xi, governed by the potential 

p 

V{Xi,...,Xp) = J2^r^l-J2^c,jkX,X^Xk + c,jkXlx]xj). (5.2) 

i= 1 ijk 

It would be interesting to investigate how this restriction affects the critical properties of 
the original model. 
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